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Abstract 

During the last few years, interest has arisen in using hght-front Tamm-Dancoff 
field theory to describe relativistic bound states for theories such as QCD. Un- 
ci^ I fortunately, difficult renormalization problems stand in the way. We introduce 
a general, non-perturbative approach to renormalization that is well suited for 
the ultraviolet and, presumably, the infrared divergences found in these sys- 
tems. We reexpress the renormalization problem in terms of a set of coupled 
inhomogeneous integral equations, the "counterterm equation." 

The solution of this equation provides a kernel for the Tamm-Dancoff in- 
tegral equations which generates states that are independent of any cutoffs. 
We also introduce a Rayleigh-Ritz approach to numerical solution of the coun- 
terterm equation. Using our approach to renormalization, we examine several 
ultraviolet divergent models. Finally, we use the Rayleigh-Ritz approach to 
find the counterterms in terms of allowed operators of a theory. 



Introduction 

First attempted in the 1950's by Tamm and Dancoff ||I|, 0, the idea of describing 
relativistic bound states in terms of a small number of particles (Fock space trunca- 
tion) ran into considerable difficulties. There were severe divergences associated with 
connected Green's functions as well as divergences associated with particle creation 
and annihilation from the vacuum and the approach was soon abandoned. In 1966, 
Weinberg noticed that, in the infinite momentum frame, creation and annihilation 
of particles from the vacuum is forbidden and the divergences associated with the 
vacuum itself are removed ||^. This is similarly true for light front quantized field 
theories. During the last few years, interest has arisen in combining these two ideas, 
Tamm-Dancoff Fock space truncation and light front quantization, in an attempt to 
describe relativistic bound states 0, H. Although divergences associated with the 
vacuum are removed, very difficult renormalization problems remain. 

In perturbation theory, divergences fall into two categories: ultraviolet (UV) diver- 
gences and infrared (IR) divergences. One generally associates UV divergences with 
particles having large energy. In the context of light-front physics, large energies arise 
from particles having large transverse or small longitudinal momenta, and thus UV 
divergences and some IR divergences have a common origin. New non-perturbative 
divergences arise because the Tamm-Dancoff truncation does not include all diagrams 
for any given order in perturbation theory. To date, these divergences in light-front 
Tamm-Dancoff calculations have only been renormalized within the context of per- 
turbation theory [Q. What we propose here is a general, nonperturbative approach 
to renormalization that naturally handles the UV and, we believe, some of the IR 
divergences found in light-front Tamm-Dancoff field theory. 

In this approach, it is shown that the renormalization problem can be rewritten 
as a set of coupled integral equations, much in the same way that the eigenvalue 
problem is written as a set of coupled integral equations. One can then solve this set 
of coupled integral equations numerically using the same techniques that are used in 
solving the eigenvalue problem itself. Although the context of light-front field theory 
is the motivation for this work, our approach will not depend on any of its special 
properties. Presumably, our approach to renormalization could be used in other, 
unrelated contexts. 

In section 1, we will introduce the renormalization problem using a simple one 



dimensional UV divergent model that has been discussed by other authors. We will 
introduce the so-called "high-low analysis" and show how it can be used to renormal- 
ize the model. In section 2 we will use the high-low analysis to examine the general 
case and derive the "counterterm equation," a set of coupled integral equations which 
relate the bare and renormalized Hamiltonians. In section 3, we will show how renor- 
malization group ideas are expressed in the context of our renormalization procedure. 
Next, in section 4, we will discuss a variational, or Rayleigh-Ritz, approach to solving 
the counterterm equation. In sections 5 and 6, we will apply our renormalization 
scheme to some simple examples. Finally, in section 7, we use the Rayleigh-Ritz ap- 
proach to find the counterterm in terms of the allowed operators of a theory. 

1. Model A 

Let us start by looking at model A, a simple toy model that has been studied by 
a number of authors P, |3, ||, H . Consider the homogeneous integral equation 



{p-E)<f){p)+g f''dp'<P{p') = (l.i: 

Jo 



with eigenvalue E and eigenvector (f){p). This is a model for a single particle of 
momentum p with Hamiltonian H{p,p') = p 6{p — p') + g. We will focus on the 
E < solution: 

, , , const ^ A , , 

Note that the eigenvalue diverges in the limit A ^ oo. Proper renormalization 
involves modifying the system to make E and (p) independent of A in the limit 
A ^ oo. Toward this end we add a counterterm Ca to the Hamiltonian. Invoking 
the high- low analysis ||T0[, we divide the interval < p < A into two subintervals: 



< p < L, a "low-momentum region," and L < p < A, a "high-momentum region," 
where the momentum scales characterized by E, L, and A are assumed to be widely 
separated. The idea is that the eigenvalue and eigenvector should be independent of 
the behavior of the system in the high-momentum region. The eigenvalue equation 
can be written as two coupled equations 



p e [0, L] ip-E)<p ip) + ig + Ca) C dp' {p') + {g + Ca) ^ dp' {p') = (1.3a) 

JO JL 



pe[L,K] [p-E)(j,{p) + {g + CK) tdp' <P{j>') + {9 + Ck) f"^ dp' (P {p') = . (1.36) 



Integrating (1.3b) in the limit L, A^ E, 



dp iP) = - . Z , /...tA L dp' ip') , (1.4) 



l + ((7 + CA)lnf io 

and substituting this expression into (1.3a), we obtain an eigenvalue equation with 
the high momentum region integrated out: 



l + ig + CA)\n^Jo 



pe[0,L] (p - E) (p) + _^^^_^Z_^ ( dp'<P{p') = 0. (1.5) 



If we demand that this expression be independent of A, 

± ( {g + CA) \ _. ._. 

rfAU + (^ + CA)lnAJ-^' ^'-^^ 

we obtain a differential equation for Ca 

dCA _ {g + Ca)' 
dk A 

Solving this equation, we are free to insert an arbitrary constant —1/A^ — In/x 



;i.7) 



At 

l-A,ln^ 



g + CA = -^ — 7^^-x- (1-^ 



Substituting this result back into (1.5), 



At I' 

1 - A„ In ^ Jo 



p e [0, L] ip-E)<P (p) + j^^^-jj I dp' ip') = (1.9) 



we see that A has been removed from the equation entirely. Using (1.8) in the original 
eigenvalue equation 



{p-E)<P (p) + ^/ ^ /^ dp' (pO = (1.10) 

I — Anln- Jo 

gives the same equation as (1.9) with L replaced by A. The eigenvalue is now, 

E = jr — — - lim E = -iie^l^^ . (1.11) 

Although the eigenvalue is still a function of the cutoff for finite A, the eigenvalue does 
become independent of the cutoff in the limit A ^ oo, and the system is properly 
renormalized. 

We can think of A^ as the renormalized coupling constant and /i as the renor- 
malization scale. In that case, the eigenvalue should depend on the choice of A^ for 
a given /x but be independent of /i itself. Suppose, for equation (1.11), we want to 
change /i to a new value, say /i'. In order that the eigenvalue remain the same, we 
must also change the coupling constant from A^ to A^i 

^e'/^^ = /i'e^/^M' . (1.12) 

In the same manner, one can write down a (3 function for A^ 

fi^A, = Al. (1.13) 



dfi 



2. The Counterterm Equation 

Using the ideas introduced in the previous section we can examine the general 
case. Throughout, we will be working with operators projected onto some Tamm- 
Dancoff subspace (finite particle number) of the full Fock space. In addition, we will 
regulate the system by demanding that each component of momentum of each particle 
lies within some finite interval. We define the "cutoff" A to be an operator which 
projects onto this subspace of finite particle number and finite momenta. Thus, for 
any operator 0,0 = AOA. Let us introduce the Hamiltonian 

H = Ho + V + Ca (2.1) 



where, in the standard momentum space basis, Hq is the diagonal part of the Hamilto- 
nian, V is the interaction term, and C\ is the counterterm which is to be determined 
and is a function of the cutoff. Each term of the Hamiltonian is Hermitian and 
compact. Schrodinger's equation can be written 

{Ho-E)<J)+{V + Ca)<J) = (} (2.2) 

with energy eigenvalue E and eigenvector 0. Our goal is to choose Ca such that E 
and are independent of A in the limit of large cutoff. 

Now we will make an important assumption: the physics that we are interested 
in, characterized by energy scale E, is independent of physics near the boundary of 
the space spanned by A. Thus, we define two projection operators, 7i and C, where 
A = H + C, TiC = CTi = 0, and H and C commute with Hq. Ti projects onto a 
"high- momentum region" which contains energy scales we do not care about, and C 
projects onto a "low-momentum region" which contains energy scales characterized 
by E. Schrodinger's equation (2.2) can be rewritten as two coupled equations: 

{Hq -E)C(j) + C{V + Ca) £0 + £ (V + Ca) 7^0 = (2.3a) 

and 

{Ho -E)n(j) + n{V + Ca) 7^0 + 7^ (V + Ca) £0 = . (2.36) 

Using equation (2.3b), we can formally solve for 7i(j) in terms of £0 

The term with the denominator is understood to be defined in terms of its series 
expansion in V. We can substitute this result back into (2.3a) 



{Ho-E)C^ + C{V + Ca)C<P + C{V + Ca) ^ {V + Ca)C<P = 0. (2.5) 



In order to properly renormalize the system, we could choose C\ such that (2.5) 
is independent of one's choice of A for a fixed C in the limit of large cutoffs. However, 
we will make a stronger demand: that (2.5) should be equal to (2.2) with the cutoff 
A replaced by C. Consider the following ansatz for Ca: 

CK = aoV + aiVFV + a2VFVFV + ... (2.6) 

where F commutes with Ti and C and the coefficients oq, ai, . . . are to be determined. 
We will discuss some specific choices for F below. Substituting (2.6) into (2.2) with 
A replaced by C everywhere and setting it equal to (2.5), order by order in V", 

(1 + ao) V + a^VCFV + a2VCFVCFV + ... 

= (1 + ao) V^ + aiVFV + a2VFVFV 

+ (1 + a,f V-,^V + ai (1 + ao) (vFV-,^V + V-,^VFv) 
Hi — riQ \ III — riQ Hi — riQ / 

+ <^ + «°''^'¥^^¥^^ + -' P.T) 

Using A = £ + 7i, we find that ag is arbitrary and the rest of the coefficients are 
uniquely determined 

a, = {-l)\l + ao)\ ^ = 1,2,... (2.8) 

provided that we can make the approximation 

V ^ V ^ VHFV . (2.9) 

E — Ho 

This is what we will call the "renormalizability condition." A system is properly 
renormalized if, as we increase the cutoffs A and C, (2.9) becomes an increasingly 
good approximation. Since the choice of ao rescales V, we absorb the factor of 1 + ao 
into the definition of V. The arbitrary scale 1 + ao corresponds to a particular solution 
of (2.7) and therefore depends on F, where F may involve an arbitrary energy scale /i. 
Thus 1 + ao is an implicit function of fi. In simple models where V is proportional to a 
coupling constant, this /i-dependent rescaling is equivalent to introducing a running 
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coupling g{fi) that runs with the arbitrary scale in F. In both models A and B 
discussed below the coupling will therefore run with fi. Summing the series (2.6), we 
can express Ca as the solution of an operator equation, the "counterterm equation," 



Ca = -VFV - VFCa (2.10) 



or, defining Va = V + Ca 



Va = V-VFVa. (2.11) 



This is our central result. In the standard momentum space basis, this becomes a set 
of coupled inhomogeneous integral equations. Such equations generally have a unique 
solution, allowing us to renormalize systems without having to resort to perturbation 
theory. This includes cases where the perturbative expansion diverges or converges 
slowly. 

There are many possible choices for F that satisfy the renormalizability condition. 
For instance, one might argue that we want F to resemble 1/ {E — Hq) as much as 
possible and choose 

^ (2.12) 



11- H( 



where the arbitrary constant /i is chosen to be reasonably close io E. In this case, 
one might be able to use a smaller cutoff in numerical calculations. 

One might argue that physics above some energy scale /i is simpler and that it is 
numerically too difficult to include the complications of the physics at energy scale 
E in the solution of the counterterm equation. Thus one could choose 

F = ^^ — ^ (2.13) 

where the arbitrary constant /i is chosen to be somewhat larger than E but smaller 
than the energy scale associated with the cutoff. The 6 function is assumed to act 
on each diagonal element in the standard momentum space basis. The difficulty with 
this renormalization scheme is that it involves three different energy scales, E, fi and 



the cutoff which might malce the numerical problem more difficult. 

3. Renormalization-Group Analysis 

We can relate our approach to conventional renormalization-group concepts. In 
renormalization-group language, Va is the bare interaction term and V is the renor- 
malized interaction term. In both of the renormalization schemes introduced above, 
we introduced an arbitrary energy scale yu; this is the renormalization scale. Now, 
physics (the energy eigenvalues and eigenvectors) should not depend on this param- 
eter or on the renormalization scheme itself, for that matter. How does one move 
from one renormalization scheme to another? Consider a particular choice of renor- 
malized interaction term V associated with a renormalization scheme which uses F 
in the counterterm equation. We can use the counterterm equation to find the bare 
coupling Va in terms of V. Now, to find the renormalized interaction term V asso- 
ciated with a different renormalization scheme using a different operator F' in the 
counterterm equation, we simply use the counterterm equation with Va as given and 
solve for V 

V = Va + VaF'V' . (3.1) 

Expanding this procedure order by order in V and summing the result, we can obtain 
an operator equation relating the two renormalized interaction terms directly: 

V' = V + V{F' -F)V' . (3.2) 

The renormalizability condition ensures that this expression will be independent of 
the cutoff in the limit of large cutoff. 

For the two particular renormalization schemes mentioned above, (2.12) and 
(2.13), we can regard the renormalized interaction term V as an implicit function 
of fi. We can see how the renormalized interaction term changes with n in case 
(2.12): 



f^^V=-V ^ V (3.3) 

rf/i {Ho - /i) 



and in the case (2.13): 



^ -^V = V5 {Ho -fx)V . (3.4) 

This is a generalization of the f3 function. 

4. The Rayleigh-Ritz Method 

It is necessary to find efficient numerical techniques for solving the counterterm 
equation. The approach we discuss here is a simple extension of a standard technique 
used in solving Fredholm-type integral equations ||lT|. Suppose we have some operator 
K; how do we measure, given V and F, whether i^ is a good approximation to V/J Let 
us define a functional J [K] such that J [K] is stationary at i^ = Va. The expression 

Tr {K - Va) {F + FVF) {K - Va) (4.1) 

is stationary at K = Va. Using the counterterm equation (2.10), (F + FVF) Va = 
FV, and (4.1) can be rewritten as 



Tr is: (F + FVF) K -Tr KFV - Tr VFK + term independent of K . (4.2) 

Thus, the functional 

J [K] = Tr K{F + FVF) K -Tr K {FV + VF) (4.3) 



is stationary at i^ = Va [0 • Let us choose a linearly independent set of Hermitian 
operators {Oj} and define 

K = J2xiOi. (4.4) 

i 

We want to find the vector x such that J [K] is stationary; this will give us the best 
approximation to Va from the subspace spanned by {Oi}. Defining the Hermitian 
matrix B and the real- valued vector c 

Bij = Tr Oi {F + FVF) O, q = Tr O, {FV + VF) , (4.5a, h) 
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the functional is 

J [TT] = x^ B X - X ■ c (4.6) 

with stationary point given by 

25 X = c . (4.7) 

Unless the number of operators is limited, this approach is not very useful. If we 
can approximate the space spanned by A using n vectors, then it will generally take 
n^ operators to span the truncated Fock space; finding x involves solving an n^ x n^ 
linear system, too large for practical applications. 

Alternatively, we can express K as an outer product of n linearly independent 
vectors {xi} 



K=Y: X.,, x^ x] ■ (4.8) 






Define the Hermitian matrices A, E^ and C: 



A, = xUf + fvf)xj 

ij Xi Xj 

a, = xHfV + VF)xj (4.9a, 6, c) 

Note that E is positive definite. The functional is now written as 

J[K] = Ti XEXA-TiCX ; (4.10) 

its stationary point is given by the matrix equation 

C = AXE + EXA. (4.11) 

One can solve for X using standard numerical techniques; the solution is basically 



equivalent to solving an n x n generalized eigenvalue problem [|14|, [15 . 
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5. Model A Revisited 

Model A is an example of what could be called a "separable system." Separable 
systems have interaction terms which can be written as an outer product of a Fock 
space vector: V = rjriK Then, the counterterm equation is 

Va = W-WFVa (5.1) 

which has the solution 

So, in this case, the counterterm equation produces an ordinary coupling constant 
renormalization. For model A in the standard momentum space basis, we have r] (p) = 
y/g where g = g (/i), and the bare interaction term is 

1 + gJo dq dq' F (g, q') 

For the choice F (g, q') = S {q — q') / {fi — q) the eigenvalue is, using a principal- value 
prescription, 

which, in the limit A ^ oo becomes independent of A, giving an expression consistent 
with our previous result, 

lim E = - 1^1 e^'^ . (5.5) 

A^oo 

We can easily verify the renormalizability condition for this model. For the Ti region, 
we restrict L < p < A. The renormalizability condition is. 
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which becomes an equahty in the hmit L, A ^ oo. Finally, we can use equation (3.3) 
to calculate the (3 function: 



which is the same result we found in (1.13). 



d 2 

lim /i 3- fi- (/i) = fi- (/i) . (5.7) 



6. Model B 

Let us consider an extension to model A that has a more complicated interaction 
term. In the momentum space basis, we define the Hamiltonian of model B to be 



^0 (p, p') =p5{p- p') V {p, p) = g 



p + p' + a 
p + p' + h 



(6.1) 



The interaction term becomes constant in the limit of large p or p', and the unrenor- 
malized eigenvalue is linear in A. 

Applying the Rayleigh-Ritz procedure, we can solve for the bare interaction term 
Va {PiP') in terms of V {p,p'). We will use (4.8) with a lattice in momentum space as 
our linearly independent set of n vectors 



Xi{p) 



1 i^A <p< ^A 
otherwise 



(6.2) 



and F {p,p') = 5 {p — p' ) / {11 — p) . Thus, the relevant matrices, A, E, and C, are: 



A 



A 5, 



^,J 



«J 



npi- fx 



A 



ViPuPj] 



nJ ij>i - fj) {Pj - fj) 



E, 



*j 



^i,j 



^* 



n 



«j 



AY 

-) yip^^PJ) 



1 



1 



Pi - /i Pj - /i 



(6.3a, b, c) 



where Pi 



i-1/2 



A. Typical numerical results for V and Va are shown in Figures 1 and 



2 forn = 100, A = 40, /x = -8.5, a = 9, 6 = 3, and ^ = -2.85. Note that Va {p,p') 
has a slightly different functional form in the regions where only one momentum is 
large and is almost zero in the region where both momenta are large. 
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Using Va, we can solve Schrodingers equation and find the renormalized eigenval- 
ues and eigenvectors. If we plot the negative eigenvalue and eigenvector as a function 
of A as shown in Figures 3 and 4, we can see that the system is properly renormalized. 

7. Local Count erterms 

An operator that is an arbitrary function of momenta in the momentum-space 
basis is generally nonlocal in the coordinate basis. In light-front physics, proper 
renormalization requires that one abandon local operators, at least for the longitudinal 
coordinates. However, when our renormalization procedure is applied to systems 
without such infrared divergences, nonlocal counterterms may still be produced. This 
is unacceptable. One possible solution to this problem is to choose only functions that 
have the desired locality when choosing a basis of linearly independent operators in 
the Rayleigh-Ritz procedure. Such a basis does not span the space of all operators, 
but the Rayleigh-Ritz procedure will, if possible, choose the coefficients so that the 
divergences are still canceled. While not conceptually pleasing, this will properly 
renormalize a system using, by construction, counterterms that are sufficiently local. 

Let us see how this works in the context of model B. Since the interaction term 
becomes constant in the limit of large p, one can renormalize model B using only a 
local counterterm. For simplicity, we will apply the Rayleigh-Ritz procedure to the 
counterterm equation for C\ (2.10) using {Oi} as our set of linearly independent 
operators where Oi is defined to be 1 in the momentum space basis. In this case, we 
want to find the stationary point of the functional 

J [K] =Tt K{F + FVF) K + Tt K {FVFV + VFVF) (7.1) 

where K = xiOi, this will give us the best possible approximation to C\ from the 
subspace spanned by {Oi}. Defining 

Ei,i = Tr Oi (F + FVF) d Ci = -Tr d {FVFV + VFVF) , (7.2a, b) 

the stationary point of (7.1) is given by 

25i,i xi = ci . (7.3) 
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In the momentum space basis, we will use F {p, p') = S {p — p') / {fi — p) and, 



Bii = -A In — + A / dp dp 



fi JfM pp' 



ci = -^J^dpdpdq — . (74^^^) 



To leading order in A, 



/ A\' 
5f A In —r - A In 



A\ A Ag{a-b)iTT^ or- /^^ + /^' 



Bii = c/A In ^ - A In ^ ^ '- <^ — + 2 Li 



2 



\ 1/^1/ l/^l 6 + 2/i [6 \ A* / J 

where Li2 (x) is the dilnarithm function, and to leading order in In ^, 

^1 = -9 (/^) + , ^^^l A • (7.6) 

in agreement with references [6] and [7]. Numerically calculating the eigenvalue as 
a function of A shows that this choice of counterterm does correctly renormalize the 
system as shown in Figures 5 and 6 where n = 100, A = 60, yU = 5, a = 8, 6 = 2.2, 
and g = —0.65. 

For a more realistic example, one would start with V written as a linear combi- 
nation of local operators and use the Rayleigh-Ritz procedure to find V\ in terms of 
those same operators. We have not yet established general conditions for the validity 
of this procedure. 

Conclusion 

We have introduced an approach to renormalization that is nonperturbative, gen- 
eral, and successfully handles the UV divergences found in light-front Tamm-Dancoff 
field theories. In this approach counterterms that remove all cutoff dependence are 
obtained from the solution of the counterterm equation. We have shown that the 
Rayleigh-Ritz procedure offers a practical way of solving the counterterm equation. 
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and we have illustrated how the renorinalization procedure works for some simple 
examples. 

In general, the Hamiltonian that one uses in field theory is not the one that 
is obtained by canonical quantization. In principle one must include all operators 
allowed by power counting. We have shown that the Rayleigh-Ritz method can be 
used to find the combination of allowed operators that provides the best solution of 
the counterterm equation. A simple example of this procedure is given. 

Although we have only examined ultraviolet divergences, the discussion that led 
to the counterterm equation should also be valid for some of the infrared divergences 
found in light-front quantized field theory. In light-front coordinates, the Hamiltonian 
for a free particle with mass m, longitudinal momentum p~^, and transverse momenta 
P± is i, i 

m'^ + pl 
12 f ^ ; • 

^ 2p+ 

Thus, if we associate the Ti region with small values of p"*", the renormalizability 
condition can be fulfilled since Hf ^ E for sufficiently small p^ . Future work will 
examine the issue of infrared divergences in detail. 

Clearly, there is a lot of work left to be done before we can successfully describe 
relativistic bound states using the light-front Tamm-Dancoff approach. However, we 
hope that the approach to renormalization outlined in this paper will provide a useful 
framework for addressing the remaining problems. 
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Figure Captions 

Figure 1. Plot of the renormalized interaction term —V {p,p') as a function of p and 
p' for model B. 

Figure 2. Plot of the bare interaction term — Va {p,p') as a function oi p and p'. 

Figure 3. Eigenvalue as a function of A along with a fit to x+y/A and the asymptotic 
value X. The fit parameters are x = —6.55 and y = 5.34. 

Figure 4. Eigenfunction p (p) as a function of p for various A. The wavefunctions 
are not normalized. 

Figure 5. Eigenvalue as a function of A along with a fit to x+y/A and the asymptotic 
value X. The fit parameters are x = —6.84 and y = 30.52. 

Figure 6. Plot of the bare interaction term — Va {p,p') as a function of p and p'. 
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